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We study the quarkonia correlators above deconfinement using the potential model with screened 
temperature-dependent potentials. We find that while the qualitative features of the spectral func¬ 
tions, such as the survival of the IS state, can be reproduced by potential models, the temperature 
dependence of the correlators disagree with the recent lattice data. 


I. INTRODUCTION 

Current and past experiments at BNL and at CERN are and have been colliding heavy ions at relativistic velocities. 
One of the major goals of these experiments is to produce a deconfined state of matter, known as quark-gluon plasma. 
In this deconfined phase the dominant degrees of freedom are the different flavors of quarks and gluons which are 
not bound inside hadrons anymore. Unlike for light quarks, due to their smaller size, bound states of heavy quarks 
could survive inside the plasma to temperatures higher than that of deconfinement, T c . In the 1980’s, however, it 
was predicted Q that cc bound states would disappear already at temperatures close to T c . The idea of Matsui and 
Satz, based on non-relativistic arguments, was that color screening in the plasma would prevent the strong binding 
of quarkonia. Therefore, the dissolution of heavy quark bound states, and thus the suppression of the J/if peak in 
the dilepton spectra, could signal deconfinement. 

It was recognized in the 1970’s )2] that what is now known as the Cornell potential (a Coulomb plus a linear part) 
provides a very good description of the quarkonia spectra at zero temperature. Much later it was understood that 
in non-relativistic systems, such as quarkonium, the hierarchy of energy scales m mv mv 2 (to and v being 

the heavy quark mass and velocity), allows the construction of a sequence of effective theories : nonrelativistic QCD 
(NRQCD) and potential NRQCD (pNRQCD) 00 . The potential model appears as the leading order approximation 
of pNRQCD gi], and thus can be derived in QCD from first principles. 

The essence of the potential model calculations in the context of deconfinement is to use a finite temperature 
extension of the zero temperature Cornell potential to understand the modifications of the quarkonia properties at 
finite temperature 0 [g]. It is not a priori clear though, whether the medium modifications of the quarkonia properties 
can be understood in terms of a temperature dependent potential. Therefore, a new way of looking at this problem has 
been developed. This new way is based on the evaluation on the lattice of the correlators and spectral functions of the 
heavy quark states 010 But why are correlators of heavy quarks of interest? As we said above, hadronic bound 
states are expected to dissolve at high temperatures. With increasing temperature such resonances become broader 
and thus very unstable and gradually it becomes meaningless to talk about them as being resonances. Accordingly, 
at and/or above T c they stop being the correct degrees of freedom. Correlation functions of hadronic currents, on the 
other hand, are meaningful above and below the transition. These can thus be used in a rather unambiguous manner 
to extract and follow the modification of the properties of quarkonia in a hot medium. 

The numerical analysis of the quarkonia correlators and spectral functions carried out on the lattice for quenched 
QCD provided unexpected, yet interesting results 010 The results suggest the following: The ground state 
charmonia, IS J/if and rj c , survive well above T c , at least up to 1.6T C . Not only do these states not melt at or close 
to T c as was expected, but lattice found little change in their properties when crossing the transition temperature. 
In particular, the masses of the IS states show almost no thermal shift ' 7,0 0 - Also, although the temperature 
dependence of the ground state quarkonium correlators was found to be small, a small difference between the behavior 
of the J/if and the r] c correlators has been identified. This difference was also not a priori expected. Furthermore, the 
lattice results indicate 0 that properties of the IP states, XcO and Xci, are significantly modified above the transition 
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temperature, and these states have dissolved already at 1.1T C . The recently available first results for the bottomonia 
states [Toj suggests, that the g b is not modified until about 2 T c , but the \bo shows dramatic changes already at 1.15T C . 

These features from the lattice calculations are in sharp contrast with earlier, potential model studies of quarkonium 
properties at finite temperature, which predicted the dissolution of charmonia states for T < 1.1T C iH More recent 
studies try to resolve this apparent contradiction by choosing a potential which can bind the heavy quark anti-quark 
pair at higher temperatures. In El the picture of strongly coupled Coulomb bound states is discussed. In 173 . fTTU *7"1~ 
the internal energy of a static quark anti-quark pair calculated on the lattice at finite temperature is used as the 
potential. It is found that the J/ip could be bound at temperatures as high as 2T c , in agreement with the lattice 
studies of the spectral functions. The quarkonium properties, however, are significantly modified in all of these 
calculations. 

In connection with the potential models two questions should be addressed: First, can the medium modifications of 
quarkonia properties, and eventually their dissolution be understood in terms of a temperature dependent potential? 
Second, what kind of screened potential should be used in the Schrodinger equation, that describes well the bound 
states at finite temperature? At present neither of these questions have been answered, since it is not clear whether 
the sequence of effective theories leading to the potential model can be derived also at non-zero temperature, where 
the additional scales of order T, gT and g 2 T are present. Therefore, in this paper we follow a more phenomenological 
approach, extending our results presented in M- We address the first question by constructing the heavy quarkonia 
correlators within the potential model and comparing them to the available lattice data. Our way of addressing 
the second question is to consider two different screened potentials and show that the qualitative conclusions do not 
depend on the specific choice of the potential. Our study of quarkonia correlators in Euclidean time is also motivated 
by the fact that although lattice calculations give very precise results about the temperature dependence of the 
quarkonia correlators (see e.g. ]§]), it is quite difficult to reconstruct the spectral functions from the correlators at 
finite temperature. 

The paper is organized as follows: Section[H]presents the model we use to construct the spectral function and analyze 
the correlators. In Section iTTTI onr results for the temperature dependence of the quarkonia properties are presented. 
Section El is devoted to the study of charmonium and bottomonium correlators in the scalar and pseudoscalar 
channels. In Section m we discuss the vector correlator which also carries information about diffusion properties. We 
summarize the results and present our conclusions in Section ED 


II. EUCLIDEAN CORRELATORS AND THE POTENTIAL MODEL 

Here we investigate the quarkonium correlators in Euclidean time at finite temperature, as these correlators are 
directly calculable in lattice QCD. To make contact with the available lattice data obtained in the quenched approx¬ 
imation, i.e. neglecting the effects of light dynamical quarks, we consider QCD with only one heavy quark flavor. 
Furthermore, we only consider the case of zero spatial momentum. 

The correlation function for a particular mesonic channel H is defined as 

Gh(t,T ) = (j H (r)j^(O)). (1) 

Here jn = q^nq, and Th = 1,7 75 , 7^75 corresponds respectively to the scalar, vector, pseudoscalar and axial vector 
channels. At zero temperature these channels correspond to the different quarkonium states shown in Table [I] Using 


TABLE I: Quarkonia channels considered. 



JH 

2S + 1 Lj 

J PC 

q = c 

-O 

II 

scalar 

qq 

3 Po 

0 ++ 

XcO 

XbO 

pseudoscalar 

q'yaq 

1 S 0 

0 ~ + 

Vc 

Vb 

vector 

Tn-.q 

3 Si 

1— 

J/1> 

T 

axial vector 

<77m75? 

3 Pi 

1++ 

Xcl 

Xbi 


the Euclidean Hamiltonian (transfer matrix) the spectral decomposition of the correlators defined in Eq. CJ at zero 
temperature can be written as fT3 | 


GH(r,T = 0) = ^|(0|iH|n)| 2 e-^ ; 


n 


( 2 ) 
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where the E n are the eigenstates of the Hamiltonian. For the few lowest lying states E n corresponds to the quarkonia 
masses, while the matrix elements |(0|jjy|n)| 2 = Ff give the decay constants. Furthermore, the Euclidean correlator 
is an analytic continuation of the real time correlator D > (t), Gh(t,T) = D^(— it, T), and the spectral function is 


cr(u,T) 


D>(u,,T)-D<(u,,T) 

27T 


1 


ImDj|(u;,r). 


( 3 ) 


Here D^(w,T), D > (ui,T) and D < (ui,T) are the Fourier transforms of the real time correlators 

( 4 ) 

( 5 ) 


D R (t,T) = m),j( 0 )]), 

D>( t ,T) = mm), o<(t,T) = mm), 


and (...) denotes the expectation value at finite temperature T. Using the above, the following spectral representation 
for Gh(t,T ) can be derived (see Appendix lAll 


G(t,T) 


duia(ijj, T) K(t, u>, T ), 


( 6 ) 


where K is the integration kernel 


T) 


cosh \u> (r - A^)] 
sinh [^] 


( 7 ) 


Recent lattice QCD calculations utilize the form © for the correlator to extract spectral functions. Results obtained 
with the Maximum Entropy Method are promising, but still somewhat controversial, and not fully reliable. We choose 
to directly analyze the correlators that are reliably calculated on the lattice. So as an input to © we need to specify 
the spectral function at finite temperature. We do this by following the form proposed in Il7l| for the zero temperature 
spectral function 


= 5Z 2M i F i5 (w 2 - M'f) + w 2 6» (w - s 0 ) /(w, s 0 ). (8) 

i 

The first term contains the pole contributions from bound states (resonances), while the second term is the perturbative 
continuum above some threshold, denoted here by sq. While we know the asymptotic behavior of the spectral function 
from perturbation theory, no reliable information is available for the explicit form of the threshold function, f(u>, so). 
Therefore, in this work we consider two simple choices for it: First, we take the most simple form of f(u>, so) = an, 
valid for free massless quarks, resulting in a sharp threshold. Then, we consider the form motivated by leading order 
perturbative calculations with massive quarks, 


/M = (a H + b H \j , (9) 

giving a smooth threshold. In leading order perturbation theory so = 2Wb jC . The coefficients (a#, bn) were calculated 
at leading order in 1181. a nd are (-1,1), (1,0), (2,1), and (-2,3), for the scalar, pseudoscalar, vector, and axial vector 
channels, respectively [32|. Comparing the results obtained with these two forms of the threshold function gives us 
an estimate on the uncertainties resulting from the simplified form of the spectral function in equation © . 

In real QCD with three light quarks it is natural to identify the continuum threshold so with the open charm 
or beauty threshold. For the case of one heavy quark only, which is considered here, the value of So is somewhat 
arbitrary. Above some energy though, the spacing between the different states is so small that they eventually form 
a continuum. Parton-hadron duality then requires, that the area under the spectral function above this energy range 
is the same as the area under the free perturbative spectral function, justifying the form given above. This also 
motivates our choice of threshold as the energy above which no individual resonances are observed experimentally. 
The parameters of the zero temperature spectral function are given in Table un 

The remaining parameters of the spectral function © can be calculated using a potential model. The bound state 
masses are given by Mj = 2m + £), where Ei are the binding energies, and m is the constituent mass of either the 
charm or the bottom quark. In order to calculate the decay constants F.;, we first write the relativistic quark fields in 
terms of their non-relativistic components. This can be done using the Foldy-Wouthysen-Tani transformation 


7 • D 

2 m 



q = exp 


( 10 ) 
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TABLE II: Parameters at T = 0. 


a 

a 

m c 

m b 

Soc 

Sob 

0.471 

0.192 GeV 2 

1.32 GeV 

4.746 GeV 

4.5 GeV 

11 GeV 


with D being the spatial covariant derivative. At leading order in the coupling and inverse mass, the decay constant 
can be related to the wave function at the origin for the S states, i.e. the pseudoscalar and vector channel EH. 

Fps = ^\Rn( 0)| 2 and F i=^\R n ( 0 )\\ (11) 

and to the derivative of the wave function in the origin for P states, i.e. the scalar and axial-vector channels 0 , 


9iV c 


Km 2 


and 


Fa = ~ 


9 N c 


K(«)l 5 


( 12 ) 


where n = 1,2,.... For the number of colors we consider N c = 3. The negative signs in Eq. o are the consequences 
of our definitions of the Euclidean Dirac matrices. For the derivation of these relations see Appendix m To obtain 
the binding energies and the wave functions we solve the Schrodinger equation 


/ 1 d 2 
\ m dr 2 


1(1 + 1 ) 

mr 2 


+ V(r) - E 


u(r) = 0, 



(13) 


with the Cornell potential 


V(r) = - h err . (14) 

r 

Here a is the coupling and a is the string tension. Their values, as well as the quark masses, were obtained in [21] by 
fitting the zero temperature quarkonium spectrum. The parameters of the zero temperature analysis are summarized 
in Table El Note that the angular momentum l in m distinguishes between the different quarkonium states shown 

in Table 0 

To model the binding and propagation of a heavy quark and antiquark in the deconfined phase, we assume that 
they interact via a temperature-dependent screened potential, and propagate freely above some threshold. Then, 
our finite temperature model spectral function has the form given in m with a now temperature-dependent decay 
constant, quarkonium mass, and threshold: 

<T(w, T) = J2 2M z (T) Fi (T) 2 S (uj 2 - M Z (T) 2 ) + m 0 u 2 6 (w - s 0 (T)) /(w, s 0 (T)) + Xs(T)q(T)u6(u ). (15) 

i 

The last term in the spectral function is important at low frequencies and contributes at nonzero temperatures. It 
is present only in the vector channel, and is due to charge fluctuations and diffusion. The function q(T) = 0 for the 
scalar and the pseudoscalar states, and q(T) = (3 T/m — 1), for the vector channel. Since we consider the case of 
zero spatial momentum, there is no distinction between the transverse and the longitudinal components in the vector 
channel. Therefore, we sum over all four Lorentz components, ay = Y2^ a v-v■ At nonzero spatial momentum, for 
u> 2 — k 2 < 0, there would also be an additional contribution to the spectral function |EJ. The charge susceptibility 
Xs in the nonrelativistic approximation has the following form 

Xs (T) = m c j^^mV 2 T 1 ' 2 e- m ' T . (16) 

The derivation of this expression, together with the full evaluation of the vector correlator at one-loop level is provided 
in the AppendixJO The relevance of this contribution in the spectral function for the possibility of obtaining transport 
coefficients from the lattice has been recently discussed in [23 ■ 

Just like at zero temperature, we calculate the binding energy and the wave function of the quarkonium states 
at finite temperature using the Schrddinger equation im For this we need to specify a temperature-dependent 
screened quark-antiquark potential. While the zero temperature potential can be calculated on the lattice, at nonzero 
temperature it is not clear how to define this quantity (for a discussion of this see 12^1 1. Some ideas on how to 
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generalize the notion of static potential to temperatures above deconfinement were presented in (20| . In the current 
analysis we use two different potentials: The screened Cornell potential above T c , first considered in Q 


V(r,T) 


_“ e -M(T) r 


H(T) 


(1 - e~^ T)r ). 


(17) 


The coupling and the string tension are specified in Table un We choose the following Ansatz for the temperature 
dependence of the screening mass fj, = [0.24 + 0.31 • (T/T c — 1)] GeV with T c = 0.270 GeV. This choice is motivated by 
the fact that we want to reproduce the basic observation from the lattice, namely that the IS state exists up to 1.6T C , 
while the IP state melts at 1.1T C . The potential is shown in the left panel of Figure ITI for different temperatures 
above deconfinement. 

In tm 12 3 the change in the internal energy induced by a static quark-antiquark pair, first calculated in )24| . 
is used as the potential. This is a debatable choice of potential, since near the transition temperature the internal 
energy shows a very large increase MEHH However, for temperatures T > 1.07T C we will also use the internal 
energy calculated on the lattice M as the potential. We parameterize this as 


V(r, T) = -- e -^ 7 > 2 + a(T)re~^ (T)r2 + C{T) (l - e~^ T)r ^ . (18) 

This parameterization is shown against the actual lattice data in the right panel of Figure E The temperature 
dependence of the parameters of ® is shown in Tables El and ED in Appendix IPl 



r [fm] 



FIG. 1: Cornell potential at T = 0 and screened Cornell potential for different T > T c (left panel); Lattice data on internal 
energy from 2 together with the fit to the data (right panel). 


In the presence of screening the potential has a finite value at infinite separation Voo(T) = lim r _ >0O F(r, T) (see 
Fig. d). The Voo(T) is the extra thermal energy of an isolated quark. The minimum energy above which the quark- 
antiquark pair can freely propagate is 2m + Voo- In what follows, we assume that this minimum energy defines 
the continuum threshold, i.e. sq(T) = 2 m + V^. Independently of the detailed form of the potential, V'oo(T) is 



FIG. 2: temperature dependence of the quark pole mass. 
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decreasing, and thus so will also decrease with temperature. This effect manifests in the temperature dependence of 
the correlators, as we will discuss in the following Section. Since at leading order the threshold of the continuum is 
2 m, we may think of sq(T)/2 as the temperature dependent pole mass, shown in Fig.|3 The decrease of Voo(T) near 
T c can thus be considered as a decrease of the pole mass. Such a decrease of the pole mass was observed in lattice 
calculations [27j, when calculating quark and gluon propagators in the Coulomb gauge. As the temperature increases 
the position of the peak will be close to the threshold, causing the phenomenon of threshold enhancement. This has 
recently been discussed for light quarks in l28j . We will not consider this in the present paper. While we assume 
that above the threshold quarks and antiquarks propagate with the temperature dependent effective mass defined 
above, quarks inside a singlet bound state will not feel the effect of the medium, and thus will have the vacuum mass. 
Therefore in the Schrodinger equation we use the zero temperature masses of the c and b quarks. Furthermore, above 
deconfinement quarkonium can also dissociate via its interaction with gluons |2?1| . This effect leads to a finite thermal 
width, which we will also neglect in the current analysis. 

Finally, in order to make a direct comparison with the lattice results we normalize the correlation function by 
the so-called reconstructed correlators 9;, [jij . This is done to eliminate the trivial temperature dependence from 
the kernel Q in the correlator ©• The reconstructed correlators are calculated using the spectral functions at a 
temperature below the critical one, here at T = 0: 

G™(t, T) = J dua(uj, T = 0 )K(t, w, T) . (19) 

The ratio G/G recon can therefore indicate modifications to the spectral function above T c . Any difference of this ratio 
from one is an indication of the medium effects. 


III. QUARKONIA PROPERTIES ABOVE DECONFINEMENT 


We first present the properties of the charmonia and bottomonia states in the deconfmed medium. These results 
are obtained by solving the Schrodinger equation m with the screened Cornell potential 113 - 




fig. 3: Temperature dependence of quarkonia masses (left panel) and amplitudes (right panel) normalized to their correspond¬ 
ing zero temperature values. 


The temperature dependencies of the masses and the amplitudes of the different quarkonia states, normalized to 
their corresponding zero temperature values, are shown in Fig. [3 The left panel of Fig. [3]shows that the bound state 
masses do not change substantially with temperature. An exception is the mass of the scalar charmonium y c o, which 
shows a significant decrease just above the transition temperature. The right panel of Fig. 0 shows the amplitudes, 
as obtained from the wave function, or its derivative at the origin. Contrary to the masses, these show a strong drop 
with increasing temperature for all the states considered. Since we neglect effects that could arise from the hyperfine 
splitting, the properties of the IS scalar rj c and pseudoscalar J/ij) states are identical. While the small shift in the 
quarkonia masses above the deconfinement temperature is consistent with lattice data, the decrease in the amplitudes 
is neither confirmed, nor ruled out by existing lattice data. 

In Fig. Q]the temperature dependence of the radii is displayed. The radii of the scalar and pseudoscalar bottomonia 
states begin to increase at much higher temperatures than their corresponding charmonia states. This is expected, 
because bottomonia is much smaller in size, and therefore survives to much higher temperatures than charmonia. As 
the temperature increases and the screening radius decreases the effective size of a bound state becomes larger and 
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T/T 

C 


FIG. 4: Temperature dependence of quarkonia radii. 


larger, as seen in Fig. 0] When the size of the bound state is several times larger than the screening radius, thermal 
fluctuations can dissociate it. This means that the spectral functions no longer correspond to well defined bound 
states, but rather to some very broad structures. Here we consider a state to be melted, when its radius becomes 
greater than 1 fm. For this reason we do not show the radius of the XcOi since it is significantly increased already at 
T c , suggesting the disappearance of this state at this temperature. Note, that the Xbo is similar in size to the r] c and 
the J/ip, and shows significant increase of its radius at about the same temperature. 

As mentioned in the previous Sections, several recent studies used the internal energy of a static quark-antiquark 
pair determined on the lattice as the potential in the Schrddinger equation to determine the properties of the quarkonia 
states at finite temperature. In the present study we also consider this possibility. The right panel in Figure E shows 
that close to T c the internal energy is significantly larger than its zero temperature value. Therefore, when used in 
the Schrodinger equation, this potential yields a large increase of the quarkonia masses and amplitudes near T c . The 
shift in the quarkonia masses is of the order of several hundred MeV, and thus disfavored by the lattice calculations 
of the quarkonia spectral functions. 

The numerical values for the quarkonia properties, together with the temperature dependence of the parameters of 
the potential are given in the Appendix |D| Tables ED and El refer to the screened Cornell potential, and Tables m 
and ED to the potential identified with the internal energy of a static quark-antiquark pair calculated on the lattice. 


IV. NUMERICAL RESULTS FOR THE SCALAR AND PSEUDOSCALAR CORRELATORS 

In this Section we discuss the results obtained for the temperature dependence of the correlators in the scalar and 
pseudoscalar quarkonia channels. We compare these to existing lattice results. 


A. Results with the Screened Cornell Potential 

We present the temperature dependence of the scalar charmonia correlator, normalized to its reconstructed corre¬ 
lator in Fig. E The left panel displays the results obtained on the lattice Q. These show a very large increase of the 
correlator at 1.1 T c , suggesting that the properties of the XcO are already modified near the transition temperature. 
The potential model calculation using (11 71) as the potential, and a sharp threshold, i.e. /(w,So) = &H in (11511 . with 
a,H = — 1 for the scalar, and 1 for the pseudoscalar, are presented in the right panel. 

The calculated XcO correlator shows an increase too, in qualitative agreement with the lattice data. Despite the fact 
that the contribution from the XcO state becomes negligible, the scalar correlator above deconfinement is enhanced 
compared to the zero temperature correlator. This enhancement is due to the thermal shift of the continuum threshold 
So(T) ' 

The pseudoscalar charmonia correlators are presented in Figure E The lattice r] c correlator (left panel) shows no 
change until about 3T c , where a decrease is detected. The potential model with a sharp threshold (right panel), 
however, yields a moderate increase in the correlator at about 0.1 fm. This feature is again attributed to the decrease 
of the continuum threshold with temperature, and is not detected on the lattice. After reaching a maximum G/G recon 
drops, due to the decrease of the amplitude Fpg(T) above deconfinement (c.f. Fig.E- 







0.4 


CD 

w 

CD 


1.8 - 


1.5 - 


1.2 


0.9 


0.75T r 

1.1T 

1.5T" 


a i 


t > * 


0.1 


0.2 

x[fm] 


0.3 



FIG. 5: Ratio of the scalar charmonia correlator to the reconstructed correlator calculated on the lattice (left panel from Q) 
and in our model with sharp continuum threshold (right panel). 
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FIG. 6: Ratio of the pseudoscalar charmonia correlator to the reconstructed correlator calculated on the lattice (left panel from 
(9|) and in our model with sharp continuum threshold (right panel). 


We now discuss the results obtained with the smooth continuum threshold ©. The temperature dependence of 
the Xco correlator obtained with the smooth threshold shows a somewhat different behavior, as seen in the left panel 
of Fig. □ After enhancement at intermediate r the correlator falls below the zero temperature one. We thus see 
that the temperature dependence of the scalar correlator is strongly affected by the continuum part of the spectral 
function. In this case too, the reduction of the continuum threshold clearly leads to the increase of the correlator at 
intermediate r. When comparing the behavior of the T] c correlator calculated with the sharp threshold (right panel 




FIG. 7: Ratio of scalar (left panel) and pseudoscalar (right panel) correlator to reconstructed correlator with smooth continuum 
threshold. 
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in Fig. 0) and with the smooth threshold © (right panel of Fig. 0 , we conclude that in the pseudoscalar channel 
although quantitatively the results are different, for the qualitative behavior only the reduction of the threshold, and 
not the exact form of the continuum matters. 




FIG. 8: Ratio of the scalar bottomonia correlator to the reconstructed correlator calculated on the lattice (left panel from 0) 
and in our model with sharp continuum threshold (right panel). 


Qualitatively similar behavior was obtained for the bottomonia states. The scalar bottomonium Xhn correlator 
shown in Fig. 0 has an increase already at 1.13T C , as determined both on the lattice (left panel from [lj|) and in 
our model calculations (right panel). Here the sharp continuum has been used. Thus, the behavior of the scalar 
bottomonium channel is very similar to that of the scalar charmonium, even though, contrary to the XcO state, the 
Xbo survives until much higher temperatures. We explain this with the fact, that the shifted continuum gives the 
dominant contribution to the scalar correlator. 

In Fig.0the behavior of the pseudoscalar bottomonium correlator is illustrated. Lattice results for this channel (left 
panel) show no deviation from one in the correlator ratio G/G recon up to high temperatures. This is considered to 
be an indication of the temperature independence of the rji, properties up to these temperatures. The potential model 
studies for the pseudoscalar correlator (right panel), yield again, a qualitatively different behavior than seen from the 
lattice. There is an increase at small r, and a drop at large r in the correlator compared to the zero temperature 
correlator. The increase is due to the threshold reduction, and the decrease is due to the reduction of the amplitude. 



FIG. 9: Ratio of the pseudoscalar bottomonia correlator to the reconstructed correlator calculated on the lattice (left panel 
from dQi]) and in our model with sharp continuum threshold (right panel). 

In summary, in this Section we have seen that G/G recon for the pseudoscalar channel is increasing above one as 
the result of the shift in the continuum threshold, and its decrease is due to the decrease in the amplitude Fps . This 
result is independent of the detailed form of the continuum. The temperature dependence of the scalar correlator is, 
on the other hand, sensitive to the continuum part of the spectral function. 
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B. Results with Higher Excited States 

In the analysis presented so far we have considered only the lowest meson states in a given channel. We have done 
so, because we wanted to make contact with the existing lattice data, in which the excited states in any given channel 
were not yet identified. In order to study possible effects due to the higher states in a given channel, we now include 
also the 2S state, and the 2S and 3S states for the pseudoscalar charmonium, and bottomonium, respectively. These 
states enter directly in the first term of the spectral function 03- Since the 2S charmonium state is melted at T c 
(see Table ITTT1 in Appendix IdI). we expect a larger drop in G/G recon compared to the previous case, where only the 
IS state is included. Similarly, the 3S bottomonium state is melted at T c , but the 2S state survives up until about 
1.75T C (see Table Hvl in Appendix iDll. The results are presented in Fig. m The left panel illustrates the temperature 
dependence of the rj c correlator obtained when both the IS and the 2S states are accounted for. As expected, we 
identify a 10%-20% reduction in the correlator, attributed to the melting of the 2S state. The lattice data on the 
pseudoscalar correlator shows no evidence for such a drop. Instead, G/G recon — 1 within errors, up to 2.25 T c . The 
right panel in Fig. m shows the rib correlator, where the melting of the 3S state near the transition, and the 2S state 
at higher temperatures, clearly induces a 10%-20% drop of the correlator compared to the case with only the IS state 
included. 



x[fm] 



FIG. 10: Ratio of the pseudoscalar correlator to the reconstructed correlator when including IS and 2S charmonium states (left 
panel) and IS, 2S and 3S bottomonium states (right panel). 


C. Results Using the Internal Energy 

We end this Section with the discussion of the quarkonia correlators obtained using as the potential the internal 
energy of a static quark-antiquark pair determined from the lattice. The temperature dependence of the correlators is 
illustrated in Fig. 1111 In the numerical analysis we used the sharp continuum. As one can see from the left panels of 
Fig-EJthe behavior of the scalar correlator is similar to the one obtained in Section flV Al Even though the \c 0 and the 
Xbo states are dissociated, above T c the correlator is always enhanced relative to the zero temperature one. As shown 
on the right panels of Fig. the enhancement in the pseudoscalar channel is larger than for the screened Cornell 
potential discussed in Section GYU As mentioned in Section EH the masses and amplitudes of the IS quarkonia 
states calculated using the internal energy as the potential, show a significant increase at T > T c . The increase in the 
amplitudes translate into a significant enhancement of the correlators with respect to their values determined with 
the zero temperature amplitudes. Besides this increase due to the amplitude, our results prove to be qualitatively 
insensitive to the detailed form of the potential, suggesting that they are based on very general physical arguments. 
Thus we tested the robustness of our results regarding the discrepancy between potential model studies and lattice 
results. 


V. NUMERICAL RESULTS FOR THE VECTOR CORRELATOR 

The vector correlator corresponds to an especially interesting channel. This is because of the fact, that the vector 
current is conserved. As mentioned in Section El (for the detailed analysis see Appendix D. this leads to an extra 
contribution in the spectral function compared to the scalar and pseudoscalar channels. This contribution (the third 
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FIG. 11: Ratio of charmonia (top panels) and bottomonia (bottom panels) correlators to the reconstructed correlators for 
scalar (left panels) and pseudoscalar (right panels) channels obtained with the potential fitted to lattice internal energy and 
sharp continuum threshold. 


term in equation m arises from the diffusion and charge fluctuations specific only for the vector state. Also, lattice 
data shows a difference in the temperature dependence of the pseudoscalar and the vector correlators. Namely, in the 
pseudoscalar channel the ratio G/G recon is equal to one up to temperatures ~ 2.25T c , and shows significant deviation 
at about 3 T c , as illustrated on the left panel of Fig. [0] In the vector channel, shown in the left panel of Fig. El 
this ratio significantly decreases from one already at 1.5T C at distances > 0.15 fm. This figure further illustrates, 
that with increasing temperature the deviation happens already at smaller distances. Since to leading order in the 
non-relativistic expansion the pseudoscalar and vector channels correspond to the same IS state, one would not expect 
the correlator of the r) c and the J/ip to behave differently. Here we conjecture, that the effects of diffusion and charge 
fluctuations make the J/ip correlator smaller than the correlator of the rj c . 


CD 


CD 


1.05 


1 


0.95 


0.9 


0.85 


0.8 



FIG. 12: Ratio of the J/ip correlator to the reconstructed correlator calculated on the lattice (left panel from f9j) and in our 
model with sharp continuum threshold (right panel). 
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The right panel of Fig. El displays the J/if) correlator from our potential model calculations. We see a rise of the 
correlator at small distances, where the continuum contribution to the spectral function is dominant, and thus again, 
the reduction of the threshold is manifest. As in the case of the pseudoscalar, the model calculations do not reproduce 
the behavior of the vector correlator obtained from the lattice. When comparing the ?7 C and the J/V' correlator from 
the model calculations, i.e. the right panels of Figs. Eland El at large distances we identify a reduction in the vector 
channel compared to the pseudoscalar channel, that results from diffusion and charge fluctuations. To better visualize 




x[fm] 


FIG. 13: Difference in the ratios of the correlator to reconstructed correlator in the vector and pseudoscalar channels for 
charmonium (left panel) and bottomonium (right panel). 


the difference between the temperature dependence of the pseudoscalar and vector correlators, in Figure El we show 
the difference between the corresponding G/G recon ratios for the charmonium (left panel) and the bottomonium (right 
panel) cases. From this figure one can clearly identify the 20% drop in the correlator relative to the r] c correlator. 
Similar to the charmonium case, for the bottomonium we also see a reduction at long distances of the T correlator 
compared to the rj & one. As expected, this reduction is more pronounced with increasing temperature, and for a given 
temperature is smaller for bottomonium than for charmonium. 


VI. CONCLUSIONS 

In this work, we presented the first detailed study of the quarkonia properties and correlators using model spectral 
functions based on the potential picture with screening, and we contrasted these with lattice data. While the potential 
model with certain screened potentials can reproduce the qualitative features of the lattice data on the quarkonia 
spectral functions, namely the survival of the IS state up to 1.6T C and dissociation of the IP charmonia state at 
T ~ 1.1T C , the temperature dependence of the corresponding meson correlators is different from that seen on the 
lattice. In general, the temperature-dependencies of the correlators show much richer structure than those calculated 
on the lattice. 

We identified several causes for the difference in the temperature dependence of the meson correlators: First, the 
properties of the different charmonia and bottomonia bound states, such as the masses and the wave functions, 
are significantly modified in the presence of screening, and will change the bound state contribution to the spectral 
function. Second, the behavior of the correlators is affected by the reduction of the continuum threshold. In a potential 
model with a screened potential the continuum threshold is related to the asymptotic value of the potential, which 
decreases with increasing temperature. Finally, the higher excited states (like the 2S, 3S, etc) are expected to melt 
at smaller temperatures above T c , inducing a large temperature dependence for the correlators. 

What are the possible implications of these findings for the picture of quarkonium binding at finite temperature? 
It is possible that the effects of the medium on quarkonia binding cannot be understood in a simple potential model. 
The effects of screening, which according to lattice calculations are present in the plasma [23|, could have time scales 
which are not small compared to the typical time scale of the heavy quark motion. In this case screening cannot affect 
significantly the properties of the quarkonia. 

On the other hand, the implications of the finite temperature lattice data on quarkonia correlators also need to 
be reexamined. It is known that the high energy part of the quarkonia spectral functions is strongly influenced by 
lattice artifacts [9], resulting in the distortion of the continuum part of the quarkonia spectral functions. Furthermore, 
the excited states were not identified in the spectral functions determined on the lattice. Based on general physical 
arguments, though, one would expect the melting of the 2S charmonia and the 3S bottomonia states close to T c , 
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together with the shift of the continuum, to induce a stronger temperature dependence of the correlators than actually 
observed in the lattice calculations. This statement is independent of whether or not the medium effects on quarkonia 
properties can be properly described in terms of potential models. In the future, therefore, this problem must be 


investigated in more detail. 
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APPENDIX A: INTEGRAL REPRESENTATION OF THE IMAGINARY TIME CORRELATOR 

Here we derive the spectral representation for the Euclidean correlator. We start from the Fourier decomposition 
of the real time two-point function 



(Al) 


The Fourier transform is given in terms of the spectral function cr(w), 

D > (u>) = 27t [1 + n(uj)} 


(A2) 


with n(u) = l/(e^ u — 1) and /3 = l/T ;3(|. The Euclidean (Matsubara) propagator is an analytic continuation of 
the real time correlator 


G(r) = D > (—ir) 


(A3) 

(A4) 



Here we insert El and obtain 



(AS) 


Now we perform a change of variables to —> — ui in the first term, and then use the identity 1 + n(u>) + n(—co) = 0 
together with the property of the spectral function, that this is an odd function in the frequency <r(—u>) = —a(to) 
(the last relation is obtained from D > {—ui) = L , which itself is derived from the periodicity in imaginary time, 
known as the Kubo-Martin-Schwinger relation We obtain 




'OO 


dbj[l + n{w)\<j{w)e TUJ 



(A6) 


Between the first and the second line above we inserted the identity 1 + n{u>) = e^“n(w). From (E® the final result 
for the Matsubara correlator is obtained: 


G(r) 



(A7) 






14 


APPENDIX B: BOUND STATE AMPLITUDES 


Here we derive the different relationships between the decay constant of a given mesonic channel and the radial 
wave function, for the pseudoscalar and vector channels, or its derivative, for the scalar and axial vector channels 
(relations (Hill and d in the main text). 

The general idea is the following: We start from the definition of the correlator for the mesonic channel H 1 


G h (t) = (0\j H {T)f H (0)\0), 


(Bl) 


where r is the Euclidean time, and the j h = qT hQ are quark currents bilinear in the heavy quark field operator, 
that distinguishes between the scalar, vector, pseudoscalar and axial vector channels via Th = 1 , 7 ^, 75 , and 7 ^ 75 , 
respectively. Using the Foldi-Wouthuysen-Tani transformation we decouple the heavy quark and antiquark fields, 
reducing in this way the relativistic theory of Dirac spinors to a non-relativistic theory of Pauli spinors. Accordingly, 


q = exp 





'jP X j + 0(l/m 2 ) 

er ■ Dvi J 


(B2) 


and 


9 = ( V’ 1 ' - X f ) - ^ ( x f cr • 5 i’ ] cr ■ D ) + 0(l/m 2 ). (B3) 

Here tl> is the quark held that annihilates a heavy quark, and x the antiquark held that creates a heavy antiquark, 

#)> = 0, <01^=0, x + |0) = 0, <0|x = 0, (B4) 

D is the covariant derivative operator, the arrows indicate its direction of action. The Euclidean 7 -matrices are given 
in the Dirac basis 


7o = 


1 0 
0 -1 


75 = 


0 1 
1 0 


7» = 


0 — ioi 

i<Ji 0 


<jj being the Pauli matrices. Note, that for simplicity we suppressed the color and spin indexes on the quark held 
operators. We can rewrite the meson correlator in the form 


Gtf(r) = <0|O*(T)Ok(0)|0>. 


(B5) 


where Oh is an operator bilinear in the quark/antiquark helds. Its particular form will be discussed in the details 
that follow. Now using 


O h (t) = e Hr O H { 0)e 


-Ht 


(B 6 ) 


and with the insertion of a complete set of states Yl n l n )( n l = where | n) is an eigenstate of the Hamiltonian, 
H\n) = E n \n), we derive a spectral representation for the correlator of the following form: 


G h {t) = Y, \(°\°H\n)\ 2 e 


-E n 


(B7) 


For large times, r —> 00 , the correlator is dominated by the state of lowest energy, and we can write 

Gh{t) ~ fne.~ ET , 


(B8) 


where fn and E are the transition amplitude and the energy of the ground state. The amplitude is related to the 
decay constant, and can be directly related to the radial wave function, or its derivative, in the origin, as it was shown 
in [ 13 . Following the steps described above let us now discuss the different channels separately. 

• We start with the pseudoscalar channel, since the calculation is simplest in this case. The current is 

jps = 9759 (B9) 

= ~ + 0(l/m) , (BIO) 
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and to leading order in the inverse mass the correlator is 

Gps(t) = (0|(xV)t(^ 1 'x)t=o|0) 

= K°lx t V’k)| 2 e“ 




(BH) 


At large r the ground state dominates, and so 

Gr, q (r ) = |( 0 | xH\Vq)\ 2 e~ En ^ 


(B12) 


where the q refers to either of the c or b quarks. The relation to the wave function in the origin is provided by 
the equation (3.12a) from 19j: 


RvA o) = 


-(0|xV|r?,) ■ 


(B13) 


Therefore, we can easily identify that 


• The scalar current is 


No 

yrj, - 27 r i^,' 


L = R n ^' 2 


(B14) 


js = qq 


= ip^ip - x*X — V ,tcr • £>X + X tcr ' Dip + x?& ■ Dip + ip^cr • £>x + 0(l/m 2 q ). (B15) 

After a straightforward, but tedious calculation, applying the relations of 03, and neglecting the disconnected 
piece, we get the following expression for the scalar correlator, to leading order in the inverse mass, 


G s (t) = - 


(2 m q ) 2 


l(°lx t<T ' Dip\n)\ 2 e~ 


■E n r 


(B16) 


Here the D denotes the difference of the derivative acting on the spinor to the right and on the spinor to the 
left, x Dip = X^ Dip — Dip- For large times the dominant contribution to the sum in mb is from the lowest 
lying state, the XgO- Utilizing equation (3.19b) from jl9) , 

\f^ R L 0 (°) = ' cr )V’IXgo), (B17) 

the relation between the amplitude and the derivative of the radial wave function at the origin is obtained for 
the XgO; with q = c,b, 

a ]\r 

' (B18) 

The m refers to the corresponding c or b quark mass. 

• Consider now the third component of the vector current: 

jv = THq 

= -W’ + <AX + X 1 Vi'0)- (B19) 

When inserting this into the correlator this yields for the diagonal components of this the following spectral 
decomposition: 

G v {t) = (0|jV( T )iv r (0)|0) 

= X!l(0|xW|n>| 2 e- B " T - (B20) 
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For r —> oo the leading contribution to this correlator comes from the vector meson V = J/ip, T, and thus 

Gj/^t) = (0|xW|Vte)>e" MvT • (B21) 

Equation (3.12b) from [p]j, with e being the polarization vector of the vector meson, is 

1 27r 


R v ( 0 )e = 

Accordingly, the transition amplitude is 


-( 0 |x T <j?/)|E(e)). 


lN r 


(0|x t (r l V’|E(e i )) = y ^i?y(0)e,, 

and summing over all polarizations e* = 3 results in 

fv = ^l-^M 0)| 2 • 

Finally, we turn to the axial vector channel. The ith component of the current is 

3 AV = 97i7s9 

= —- X t CT i X 

i 


(B22) 


(B23) 


(B24) 


+ 


2 m 0 


V’ T o' J cr • £)y - x^er • D4> + X f(T ' D(Ji4> - ’ D<JiX + 0(1 /ml) ■ 


(B25) 


When evaluating the correlator we again make use of the cancellations of CEO), and the anticommutativity of 
the quark field operators. Furthermore, we also apply the relation between the non-commuting Pauli matrices, 
(TjOj = Sij + icijk&k, and the tensor notation of a vector product, (a x b)i = £jThe resulting form of the 
correlator, when considering only Xqi, with q = c,b, and neglecting the disconnected part, the grounds state, is 


GW(r) = -- 


X9lV '' ( 2 m g ) 2 


|( 0 |x t (-i£) x er)V>|x g i(e))| 2 e Ex “ lT 


(B26) 


Equation (3.19c) from E3 relates the transition amplitude to the derivative of the wave function at the origin: 

(B27) 


Q AT 1 

(°) e = -yf (°lx + (— X o-)V’lx«i(e)) ■ 

We identify then the decay constant of an unpolarized axial vector channel to be 


fx ql = -^K al m 2 - 


(B28) 


APPENDIX C: VECTOR CORRELATOR AT ONE-LOOP LEVEL 


Here we derive the meson spectral function as obtained from the imaginary part of the retarded correlator. We 
will focus on the vector channel, since the spatial and the temporal component of its spectral function contain terms 
relevant for the analysis of the transport properties. The spectral function in the scalar, pseudoscalar and axial 
channels can be derived in an analogous manner. 

The correlation function of a meson can be evaluated in momentum space, using the quark propagator and its 
spectral representation. Accordingly, 


Xh(u,p) = N C T ^2 

n 


d 3 k 

(M 


rTr 


T H S F (k 0l k)r2sUk' 0 ,k') 


(Cl) 


where N c is the number of colors, and H refers to a given mesonic channel, defined by the operator T# = 1, 7 ^, 75 , 7^75 
for the scalar, pseudoscalar, vector, and axial-vector channels, respectively. Here k' 0 = ko — u>, k' = k — p and 
ko = iu n = i(2n + l)irT. The quark propagator has the following form -3Qj 

pi/T poo 

Sp(ko, k) = — ( 70 /co — 7 • k + m) / dre k ° T / dupFipj) [1 — n(w)] e _a,T . (C2) 

J 0 J — 00 
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Here pf(lo) = 5{to 2 — w%) is the free quark spectral function, with to 2 = k 2 + m 2 for a quark of mass to. After inserting 
El into EB and evaluating the integrals over the imaginary time r, one obtains 


Xh{u,p) = N c T^2 

n ' 

= N ^Y 


d 3 k 1 


1 


( 2 tt ) 3 u> 1 - kl to% - k’ 2 

d 3 k 1 1 

( 2 tt ) 3 to 2 - k 2 to 2 , - k’ 2 


Tr 


+ m)T ] H {Yk ,v + to) 

4 [kPh"' + k u k'» + g^(m 2 - k ■ k')\ 


where for the second line we specified the vector correlator, Tjj = 7 ^. 

First, let us evaluate the spatial part of El, i-e consider p, v = i,j and sum over the components: 


V- ^ / d 3 k 1 1 

Y^XiMP) = natYJ ^ fLol _ k 2^_ k ,2, 


l=J 


Y | kik'j + kjkl — Sij (m 2 — (kok' 0 — k ■ k’ 


i=j 


= natY 


d 3 k 1 


(2 n) 3 co 2 -k 2 to 2 ,-k' 2 


[2 k 2 + pk cos 6 + 3(wfc — Aq) — 3A: 0 w] . 


(C3) 


(C4) 


Here we inserted p • k = pk cos 6. Using standard techniques [2fl| the Matsubara sums can be evaluated: 


T E ( „ 2_;.2 - n 2 n * 0 > 


~ k o 


T E 


t E 


1 


2 w k 
1 


K-*§)(“&-fc& 2 ) 


{ul-k 2 ){u>l,-k' 2 ) 


4,LOkLOk' 


(1 - n k - n k ') - 


1 


1 


1 

4 LOk 


(1 - n k - n k ') - 


LO - LOk - LOk' LO + Wfc + LOk' 


1 


( 1 

1 

Vw + LOk - LOk' 

LO - LO k + LO k ' J _ 


1 


+(rifc - rife/) ( — 


_ + _ 
tO — Wfc — LOk' LO + Wfc + LOk' 

1 1 

+ 


W + LOk - LOk' to - LOk + LOk' 


(C5) 


where rik = l/(exp(w/ c /T) + l) is the Fermi distribution function, with tok = Vk 2 + to 2 . In the following we introduce 
a small momentum approximations p -C fc, according to which 


dn k kp 

rik'-rik — -- -cost*, 

ClUJk LUk 

nk< ”h nk — ‘Znk > 

pk a 

LOk' - LOk — -COS 0 , 

LOk 

LOk' + LOk — 2 LOk ■ 


(C 6 ) 

(C7) 

(C 8 ) 

(C9) 


The spectral function is given by the imaginary part of the momentum space retarded correlation function, \^{to,p) = 
Xij(u>+ ie,p), 


Y (j h(u,p) = — Im E ( w ’ p) ■ 


(CIO) 




Im- 


10 — E + ie 


We apply the relation 


= —ttS(lo — E ), 


(Cll) 

(C12) 
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and obtain 

^2<Ju(u,p ~ 0) = N c 


d 3 k 

( 2 ^F 


\( 2k- + pk cos 9 3uA 

s - 2 -1-(1 - 2 n k )S(co - 2us k ) 

IV u k J 


= Nr. 


d 3 k 


2k 2 +pk cos 9 dn pk 
lo\ duj k uj k 
3 u) dn pk 
iOk dui k UJ k 

2k 2 + pk cos 9 3 Q 


cos 9 


pk 


Wfe 


S ( u) -cos 9) — 5 ( uj H-cos 


pk 


Wfc 


COS0 


pk 


Wfe 


S ( u) -cos 9 \ + 5 [w + — cos 9 


pk 




2k 2 dn pk 


H-) (1 - 2n k )S(u - 2uj k ) 

u>k. 


pk 


—2— 7T-, -cos 96 [ lo -cos0 


ujf, duj k 0Jk 


Uk 


(C13) 


(C14) 


In the integrand x5(lo — x) —> ujS(lo) for the limit of zero momentum, x —> 0. The result for the spectral function at 
zero momentum is then 


'£2<T ii (u),p = 0) = 


N, 


c ,-,2 


8 tt 2 




4 TO 2 


tanh 


u) 


4 m* 


- 2 N c u5(Cd) 


d 3 k dn 2k 2 


(27r) 3 dwfe w 


(C15) 


The first term is the leading order perturbative result for the continuum, also obtained in m- The second term 
provides a temperature-dependent contribution to the correlator. We evaluate this term in a non-relativistic approx¬ 
imation at the end of this Appendix. 

Next, we evaluate the temporal part of the correlator. Inserting p = is = 0 in o results in 


Xoo (w,p) = N c 4T^2 


d 3 k 


1 


1 


( 2 tt ) 3 u 2 - k 2 uj 2 , - k' 2 L 

d 3 k_ 

3 v fc ~ a 0 w k> ~ ^0 

The corresponding component of the spectral function at small momentum is 


2kok' 0 + (jn 2 — (k 0 k' 0 — k ■ k')^ 


r (pi~ i i 

= N ‘ iT E / (2j) 3 _ tg ^ - fcg. [ 2 “* ~ * k ~ (“»- fc o 2 ) - M 


1 . 


<T 00 (u;,p ~ 0) =-Imxoo(w,p ~ 0) 


= AV 


d 3 fc 

( 2 ^F 


/ 2ul-pkcos9 w\,. „ .,, „ . 

---o- 1 - (!- 2n fc )5(w - 2 w fc ) 

V UkJ 

dn pk (_ pk 

+4--cos 05 ( a;-cos ( 


LOk 

In the limit of zero momentum the contribution from the first term vanishes, and the final result is 

d 3 k dn k 


cr O o(w,p = 0) = 4N c lu5{u>) j 


(27r) 3 dujk 




(C16) 

(C17) 

(C18) 

(C19) 

(C20) 


Here we identified the static charge susceptibility, Xs ? which can be evaluated explicitly. For heavy quarks the Fermi 
distribution function is well approximated by the Boltzmann distribution, n k = exp (— u> k /T ), and thus 


Xs = -4 N c 


= 4 Nr 


d 3 k dn k 
(27t) 3 du>k 


2ir 2 T 


dkk e T 


= 4 N r . 


1 


2tt 2 T, 
1 
27T 2 


duju>(i 


(J 2 - 


= 4 N c ^m 2 K 2 


(C21) 

(C22) 
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where K 2 is the modified Bessel function of the second kind. Since m/T 1 the following approximation can be 
applied 


K v {z) 



for z > 1, 


yielding for the static susceptibility 


x - = 4JVi (syfl m3,2rI/2 '" m/T - 


(C23) 


(C24) 


Finally, we now evaluate the second contribution to cm while incorporating a non-relativistic approximation for 
the heavy quarks. This means n k = exp (— uik/T) and k/m <C 1. Thus 


^2(ni(u),p = 0 ) 
i 


—4:N c ijj5(u)) 


4N c lu5(lo) 


d 1 2 3 4 k dn k 2 
(27r) 3 du k 


2tt 2 T 


o-m/T 


dk 


0 -k 2 /(2mT) 


k 2 


4 Nr. 


2tt 2 Tto 2 
3QS(u>) ( 4TV, 


p OO 

u5{u>)e- m/T / dkk A e~ k2/(2mT) 

J 0 


( 2 t r ) 3 / 2 


m 3/ 2r l/2 e - m /T 


T 

m 


T _ 

3 x s —u5(lu) . 
m 


(C25) 


APPENDIX D: NUMERICAL RESULTS AS TABLES 


Here we provide the results of the numerical calculations for the properties of the charmonia states, in Table um 
and bottomonia states, in Table hyi together with the temperature dependence of the screening mass in the screened 
Cornell potential, as well as the parameters of the potential fitted to the lattice internal energy and the charmonia 
properties in Tables El and eh 


TABLE III: Results for the cc mesons. 


T/T c 

M[GeV] 

state 

M[GeV] 

<r 2 > 1/2 [fm] 

|P(0)| 2 (S); |R'(0)| 2 (P) 

0 

0 

IS 

3.07 

0.453 

0.735 GeV 14 



2S 

3.67 

0.875 

0.534 GeV 3 



IP 

3.5 

0.695 

0.06 GeV 5 

1.1 

0.25 

IS 

3.044 

0.574 

0.468 GeV 3 



2S 

3.338 

1.591 

0.137 GeV 3 



IP 

3.294 

1.124 

0.012 GeV 5 

1.25 

0.314 

IS 

3.031 

0.633 

0.395 GeV 3 



IP 

3.23 

1.58 

0.005 GeV 5 

1.5 

0.395 

IS 

3.012 

0.73 

0.311 GeV 3 



IP 

3.132 

1.776 

0.001 GeV 5 

1.75 

0.472 

IS 

2.992 

0.89 

0.23 GeV 3 

2 

0.55 

IS 

2.968 

1.2 

0.15 GeV 3 


[1] T. Matsui and H. Satz, Pliys. Lett. B 178, 416 (1986). 

[2] E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane and T. M. Yan, Phys. Rev. D 21, 203 (1980). 

[3] G. S. Bali, Phys. Rept. 343, 1 (2001) arXiv:hep-ph/0001312 . 

[4] N. Brambilla et al, Quarkonium Physics , CERN Yellow Report, hep-ph/0412158 

[5] F. Karsch, M. T. Mehr and H. Satz, Z. Phys. C 37, 617 (1988). 



























20 


TABLE IV: Results for the bb mesons. 


T/Tc 

M[GeV] 

state 

M[GeV] 

(r 2 > 1/2 [fm] 

|i?(0)| 2 (S); |R'(0)| 2 (P) 

0 

0 

IS 

9.445 

0.225 

9.422 GeV 3 



2S 

10.004 

0.508 

4.377 GeV 3 



IP 

9.897 

0.407 

1.39 GeV 5 



3S 

10.355 

0.74 

3.42 GeV 3 



2P 

10.259 

0.65 

1.74 GeV 5 

1.1 

0.25 

IS 

9.52 

0.243 

8.346 GeV 3 



2S 

9.945 

0.641 

2.632 GeV 3 



IP 

9.893 

0.501 

0.738 GeV 5 



2P 

10.092 

0.95 

0.54 GeV 5 

1.25 

0.314 

IS 

9.536 

0.249 

8.07 GeV 3 



2S 

9.924 

0.695 

2.195 GeV 3 



IP 

9.885 

0.541 

0.574 GeV 5 



2P 

10.041 

1.13 

0.32 GeV 5 

1.5 

0.395 

IS 

9.553 

0.259 

7.644 GeV 3 



2S 

9.895 

0.806 

1.676 GeV 3 



IP 

9.872 

0.603 

0.404 GeV 5 

1.75 

0.472 

IS 

9.569 

0.267 

7.257 GeV 3 
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TABLE V: Results for the cc mesons using the potential fitted to lattice internal 
energy. 
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|i?(0)| 2 (S); |R'(0)| 2 (P) 

0 

0.18 

0 

0 

IS 

3.177 

0.500 

0.441 GeV 3 





2S 

3.733 

0.923 

0.370 GeV 3 





IP 

3.534 

0.733 

0.034 GeV 5 

1.07 

0.0417 

0.123 

1.467 

IS 

3.522 

0.399 

0.863 GeV 3 





IP 

4.059 

0.791 

0.072 GeV 5 

1.13 

0.0484 

0.155 

1.09 

IS 

3.448 

0.458 

0.688 GeV 3 

1.25 

0.074 

0.193 

0.747 

IS 

3.448 

0.636 

0.411 GeV 3 


[6] S. Digal, P. Petreczky and H. Satz, Phys. Rev. D 64, 094015 (2001) arXiv:hep-ph/0106017 . 

[7] T. Umeda, K. Nomura and H. Matsufuru, Eur. Phys. J. C 39S1, 9 (2005) arXiv:hep-lat/0211003 . 

[8] M. Asakawa and T. Hatsuda, Phys. Rev. Lett. 92, 012001 (2004) arXiv:hep-lat/0308034 . 

[9] S. Datta, F. Karsch, P. Petreczky and I. Wetzorke, Phys. Rev. D 69, 094507 (2004) arXiv:hep-lat/0312037 . 

[10] K. Petrov, A. Jakovac, P. Petreczky and A. Velytsky, arXiv:hep-lat/0509138 

[11] E. V. Shuryak and I. Zahed, Phys. Rev. C 70, 021901 (2004) arXiv:hep-ph/0307267 ; Phys. Rev. D 70, 054507 (2004) 
arXiv:hep-ph/0403127 . 

[12] C. Y. Wong, arXiv:hep-ph/0408020 

[13] W. M. Alberico, A. Beraudo, A. De Pace and A. Molinari, arXiv:hep-ph/0507084 

[14] M. Mannarelli and R. Rapp, arXiv:hep-ph/0509310 

[15] A. Mocsy and P. Petreczky, Eur. Phys. J. C 43, 77 (2005) arXiv:hep-ph/0411262 ; arXiv:hep-ph/0510135 

[16] I. Montvay, G. Munster, Quantum Fields on a Lattice, Cambridge University Press, 1996. 













































21 


TABLE VI: Results for the bb mesons using the potential fitted to lattice internal 
energy. 
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